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We show how to achieve optical, spin-selective transitions from the ground state to excited orbital 
states of group-V donors (P, As, Sb, Bi) in silicon. We consider two approaches based on either 
resonant, far-infrared (IR) transitions of the neutral donor or resonant, near-IR excitonic transitions. 

For far-IR light, we calculate the dipole matrix elements between the valley-orbit and spin-orbit 
split states for all the goup-V donors using effective mass theory. We then calculate the maximum 
rate and amount of electron-nuclear spin-polarization achievable through optical pumping with 
circularly polarized light. We find this approach is most promising for Bi donors due to their large 
spin-orbit and valley-orbit interactions. Using near-IR light, spin-selective excitation is possible for 
all the donors by driving a two-photon A-transition from the ground state to higher orbitals with 
even parity. We show that externally applied electric fields or strain allow similar, spin-selective 
A-transition to odd-parity excited states. We anticipate these results will be useful for future 
spectroscopic investigations of donors, quantum control and state preparation of donor spin qubits, 
and for developing a coherent interface between donor spin qubits and single photons. 

PACS numbers: 71.55.Ak, 78.55.Ap, 03.67.-a, 71.15.-m 


I. INTRODUCTION 


The optical spectroscopy and control of shallow group- 
V donors (P, As, Sb, Bi) in silicon has a long history®® 
Recent work has focused on achieving optical control 
of the group-V donors for applications to quantum 
information 3 and far-infrared (IR) lasers.- One common 
approach takes advantage of the excited orbital states of 
the neutral donors (D° states), which have hydrogen-like 
s-, p-, d-, etc. orbitals.' 5 ^ Resonant excitation into these 
states allows one to turn on and off interactions between 
donors by exciting and de-exciting the system®® or 
create population inversion for lasingP® 

The second class of approaches makes use of the neu¬ 
tral, donor bound exciton state consisting of a hole bound 
to the doubly occupied donor (D°X states )P This state 
is appealing to use for optical control because it has a 
long lifetime, a large spin-orbit splitting, and an eas¬ 
ily accessible transition frequency in the near-infrared.^ 
It has been successfully employed for electron and nu¬ 
clear spin state preparation and measurement and spin- 
to-charge conversion for spin readout of P and Bi donor 
ensembles JE303] 

In this Article, we describe several approaches to 
achieve optical, spin-selective excitation from the ground 
state to the excited orbital states. Figure [l] illustrates the 
energy scales and terms relevant for optical spectroscopy 
and control of the group-V donors. We first analyze 
control schemes based on near resonant, far-IR transi¬ 
tions. These transitions are most promising for Bi donors 
because of their large spin-orbit and valley-orbit inter¬ 
actions. Using multi-valley effective mass theoryp®®^ 
we calculate the dipole matrix elements between the 
spin-valley-orbital states for all the donors and estimate 
the maximum rate and amount of electron-nuclear spin- 



FIG. 1: Energy scales in optical spectroscopy of group-V 
donors in silicon. Here Ai, T 2 and E refer to irreducible 
representations of the tetrahedral group, J is the total angu¬ 
lar momentum quantum number, F is the hyperfine quantum 
number, and I is the total nuclear spin quantum number. 


polarization achievable with this technique. We then con¬ 
sider optical approaches in the near-IR using the D°X 
state. In this case, spin-selective excitation is possible 
for all the donors by exploiting a two-photon A-transition 
from the ground state to higher orbitals with even par¬ 
ity (e.g., 2s, 3d). Using externally applied electric fields 
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or strain, we show it is also possible to achieve similar, 
spin-selective A-transitions to odd-parity (e.g., 2 p) ex¬ 
cited states. 

The paper is organized as follows. In Sec. [IT] we review 
the effective mass theory (EMT) used to describe the 
donor orbital wavefunctions and the donor bound exciton 
state. In Sec. m we use EMT to calculate parity forbid¬ 
den dipole matrix elements between spin-valley-orbital 
states of the donors. To take into account the large cor¬ 
rections from the donor impurity potential for the Is 
states, we use the recently developed variational solu¬ 
tions provided in Ref. HD Using these calculated dipole 
elements, we consider the achievable amount and rate 
of electron-nuclear spin polarization via optical pump¬ 
ing through the spin-orbit split 1 sT 2 states. In Sec. |TV| 
we consider several different schemes to achieve spin- 
selective excitation to the excited orbital states of the 
donor through the D°X state. We calculate the relative 
strength of the dipole transition from the D°X state to 
the even-parity states by using a Hartree-Fock variational 
solution to the D°X ground stated We then show how 
to excite the odd-parity, 2 p states through two-photon 
transitions using a static electric field or externally ap¬ 
plied strain. Throughout this work we use simple sym¬ 
metry arguments and variational solutions to obtain or¬ 
der of magnitude estimates for these effects. A detailed 
and accurate understanding of these effects will require a 
combination of more refined measurements and ab-initio 
theory. 


II. EFFECTIVE MASS THEORY 


In this section we outline the formulation of the effec¬ 
tive mass theory and variational solutions to the donor 
wavefunction orbitals. The donor atom is a point-like de¬ 
fect in the lattice, while the bound electron wavefunction 
extends across hundreds of lattice sites. As a result, the 
combined system probes the silicon lattice at both long 
and short wavevectors. Early on it was realized that 
the excited orbital states cam be understood within a 
hydrogen-like model and have a universal spectrum.® The 
electron in the ground state, however, has a large overlap 
with the donor atom, resulting in a large, donor-specific 
valley splitting, i.e., chemical shift® The interaction with 
this core potential results in sufficient complexity that a 
first-principles understanding of these effects is only now 
being developed through a combination of scanning tun¬ 
neling microscopy measurem ents, atomistic simulations, 
and variational methods! 1 8 I 22 I 23 I 25 IUH1 

In what follows, we assume the validity of electron-like 
quasiparticle excitations on the semiconductor vacuum 
with a local crystal potential and a static, isotropically 
screened Coulomb interaction. This approximation has 
been very successful in semiconductors and results in an 
effective Schrodinger equation for the donor electron in 


the presence of an impurity of the form 


E 


ft 2 V 2 

2toq 


+ V°{r) + U(r) 




(1) 


where H is Planck’s constant, m 0 is the electron mass, 
V°(r) is the periodic potential of the undoped Si lattice, 
U (r) is the attractive potential of the impurity, and E is 
the electron energy. 

Si is an indirect band gap semiconductor and, con¬ 
sistent with the tetrahedral symmetry, the conduction 
band has six degenerate valleys located along the crys¬ 
tallographic (100) directions at the points fc^o = koi with 
% = ±x,±y, ±z and k 0 « 0.85 x 2tt/ a S i ( a S i = 0.543 nrn 
is the lattice spacing). 

The eigenstates of Eq. (|T]) take the form 

W = [ ^kMk^ + k^o), ( 2 ) 

(r\k^ + k^o) = w 0 (*U + r)e l(fc '‘ +fc ' i ° ) ' r , (3) 


where /i = ±cc, ±y, ±z sums over the conduction band 
minima and Uq is the Bloch wavefunction, which is in¬ 
variant under lattice translations. The effective mass ap¬ 
proximation assumes <C k ^o so that we can replace 
+ k fl0l r) « rio(fc M o,r). In this case 


(r\ip) « ^ a^F t _ l (r)u^(r)e lkfJ - 0 " r , (4) 


where F fl is the envelope function. Far away from the 
impurity U(r) reduces to the Coulomb potential and F fl 
satisfies the effective mass equation 


EF^r) 
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( 5 ) 


where 9 2 j_ = <9 2 + <9 2 (similarly for other y), « 

0.916(0.191) mo is the effective mass in the direction par- 
allel(perpendicular) to /i, mo is the electron mass, eo 
is the dielectric constant, e is the electron charge, and 
e ~ 12 is the static dielectric constant of Si. 

Equation <[5j) suggests that the ground state of the 
donor is six-fold degenerate; however, this degeneracy 
is broken by intervalley coupling induced by the donor 
atom. These effects can be self-consistentl y incl uded in 
the EMT through the multi-valley equation 15 16 


EF^ = +Y,K u " ei(k “ 0 ~ k * 0) ' rU ( r ) F "’ 


U(r) = - 
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+ U cc (r ) 


( 6 ) 

( 7 ) 


where is the anisotropic kinetic energy operator from 
Eq. ^ and U cc is the so-called “central cell” potential, 
which takes into account deviations from the Coulorn- 
bic potential in the vicinity of the donor. Several ap¬ 
proximation schemes have been developed to extract U cc 
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based on fitting the potential by comparing experimen¬ 
tally measured quantities such as energy splittings or 
the hyperfine coupling to the same quantities extracted 
from variational solutions to Eq. ©>• The solutions to 
Eq. @ will only be reliable if the resulting F /1 (k) re¬ 
main strongly localized around k^Q. Provided this con¬ 
straint is satisfied, the multi-valley EMT is a power¬ 
ful computatonal approach which has provided insight 
into the electronic structure and relaxation rates of the 
donors hyperfine and quadrupolar i nterac tions of 
the donor nucleus E2JI2I static St ark e ffects] 20 1 21 1 and ex¬ 
change coupling between donors p2M! 

Although Eq. breaks the valley degeneracy, the ef¬ 
fective Hamiltonian still commutes with the tetrahedral 
symmetry group T&- As a result the eigenstates of Eq. 
([6]) split into irreducible representations of Td . The char¬ 
acter group for the six valley sites in the reciprocal lattice 
splits into Xvaiiey = Ai + E + T 2 . The lowest energy state 
of Eq. ([ 5 ]) is a ls-like state. The Kohn-Luttinger varia¬ 
tional solution to F z away from the central cell is given 

b 


p z ~ 1 p -- v /(a 2 +y 2 )/a 2 +z 2 /b 2 (g) 

\/ nab 2 

and similarly for the other /i, where a and b are varia¬ 
tional parameters. The six Is states have the same sym¬ 
metry as Xvaiiey and, in the orbital space, the represen¬ 
tations are given by 

OLA\ ^/^(l) 5 

(*E = {*( 1 , 1 ,—1,-1,0,0), ^=(1,1,1,1,-2,-2)J, 
a T2 = |-^(1,-1,0,0,0,0), -^(0,0,1,-1,0,0), 
-^=(0,0,0,0,1,-1) | 

where the six entries correspond to the 
(+X, —X, +Y, — Y, +Z, —Z) valleys, respectively, and 
(r| ipn) = Y^i_i a nFn, f i u fi elk ‘*°" r ■ These irreducible repre¬ 
sentations are also important for the 2s states, which 
have a large overlap with the donor nucleus; however, 
for states such as 2 p 7 which vanish at the donor site, the 
valleys remain decoupled and, to a good approximation, 
F njI is described by Eq. ^ . 


III. FAR-IR SPIN-VALLEY-ORBITAL 
CONTROL 

In this section we show how to achieve spin-selective 
excitation of the donor ground state by exciting the sys¬ 
tem with far-IR light through the spin-orbit split 1 sT 2 
states. We calculate the dipole moment for this transition 
and estimate the achievable electron spin polarization. 
This approach should be most effective for Bi donors due 



P 

As 

Sb 

Bi 

Ei S T 2 -r 7 ~ E lsAl (meV)“ 

11.7 

21.1 

9.7 

38.1 

EisT 2 ;r s — E 1 s a 1 (meV)“ 

- 

- 

10.0 

39.1 

E 2 p 0 — Ei sAl (meV)° 

34.1 

42.3 

31.2 

59.5 

Eip± — E 2 P0 (meV)“ 

5.1 

5.1 

5.1 

5.1 

A (meV) 

Theory 

0.02 6 

0.09^ 

0.34 6 

1.03 6 

Exp. 

- 

- 

0.29“ 

1.0“ 

fhs Al ,lsT 2 (D) 

Theory 

0.02“ 

0.04“ 

0.02“ 

0.05“ 

Exp. 

- 

- 

- 

~ l d 

fils,2p 0 (D) 

Theory 

31“ 

25“ 

34“ 

15“ 

Exp. 

13 f 

10-30 9 

- 

- 

f-iis,2p± (D) 

Theory 

59“ 

45“ 

65“ 

28“ 

Exp. 

49 h 

•C 

CO 

CO 

- 

- 


“Exp., Ref. El 
^Theory, Ref. 1301 

“Theory, this work. 

“Exp., Ref. PTT1 .se e Ref. |29| pg. 180 for reproduction of data. 
“Theory, Ref. 1351 
^Exp., Ref. [8] 

9 Exp., Ref. 1331 
h Exp., Ref. El 

TABLE I: Table of parameters relevant for far-IR spin-valley- 
orbital control. Here y? ab = | JA =a , . |(5|i|a}| 2 is the average 
dipole moment given in Debye. 


to the large nuclear mass and, thus, spin-orbit splitting 
induced by Bi. 

The spin-orbit splitting of the excited states of the 
donor is negligible due to the strong dielectric screening 
and weak spin-orbit coupling in silicon. However, in the 
vicinity of the donor nucleus, the Is states have a strong 
interaction with the donor nucleus, which can enhance 
the spin-orbit splitting. In particular, 


IF 


h 2 

2toqC 2 


VC/ x p ■ s 


h 2 1 dU cc 
2toqC 2 r dr 


(9) 


where p, s and L are the electron momentum, spin, and 
angular momentum, respectively. The spin-orbit cou¬ 
pling matrix elements in the 1sT 2 manifold are given by— 


A =2^< T -I 

3 h 2 


( 10 ) 


2 TOgC 2 


d 3 rF*F r u*u x - 


1 dlL 


dr 


■hoy sin koy cos k$x 


The 1 5 X 2 states have a spin-one representation with re¬ 
spect to L and split into effective spin-1/2 and spin-3/2 
manifold separated by A. The eigenstates | J,mj) are 
characterized by the usual total angular momentum J 
and and ^-angular momentum mj such thaP^ 


J= 1/2 : 11/2,1/2) = y/2j3\T+ i) - v^l^o t>, (H) 

J = 3/2 : 13/2, 3/2) = \T + t), (12) 

|3/2,1/2) = y/l/3\T+ i) + v^lTo f), (13) 
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( a )j = 3/2 mj = 3/2 O 3 ) ^ (c) <J + 



F = 1-1/2 


FIG. 2: Selection rules and Clebsch-Gordon coefficients for (a) left-hand-circularly polarized light er_, linearly polarized light 
along the 2 -axis tt z , and (c) right-hand-circularly polarized light a+. We also show the two hyperfine manifolds, however, this 
energy splitting is typically smaller than the IsTi linewidth so is not resolvable in these transitions. 


Here | T m ) are the projections of L z eigenstates into the 
1sT 2 space. 

Taking F ^ to be of the form of Eq. (l8| with a and b 
given by the Kohn-Luttinger values^EH gives an estimate 
that the spin orbit interaction should be reduced from 
the atomic value by a factor of 10 -3 — 10” 4 [.‘14j . This 
rough estimate is consistent with the calculated values 
(P, As, Sb, Bi) and the measured values (Sb, Bi) shown 
in Table |U 


A. Parity-Forbidden Transitions 


Taking advantage of these spin-orbit split states for 
quantum control requires optical addressing of these 
transitions. Unfortunately, the IsAi —> 1 sT 2 transition 
is neither dipole nor Raman allowed within the single val¬ 
ley EMT. This is because the Is states are parity eigen¬ 
states within the EMT with eigenvalue +1, while dipole 
transitions should change the parity when it is a good 
quantum number. At the same time, the absence of Ra¬ 
man transitions follows because the 1sT 2 states are anti¬ 
symmetric combinations of opposing valley states, while 
the IsAi states are symmetric combinations; thus a very 
large wavevector ~ ko is required to induce a transition 
between these states. 

Although parity is preserved within the single valley 
EMT, parity is not a good quantum number within the 
tetrahedral group, thus, within the multi-valley EMT 
these transitions become dipole allowed. The dipole ma¬ 
trix element is given by 


li = (T 0 |2|1sAi) 


i 

71 


d 3 r(F^ 2Z F AuZ zsm2k 0 z 


+ AF^ 2<z F AliX z sin k 0 z cos k 0 x), (14) 


where we have approximated the product of Bloch func¬ 
tions = 1 for all v and fi. Similarly (T m |a; ± 

iy\lsA\)/\/2 = fiS rn ±. The polarization selection rules 
for the IsAi to 1 sT 2 triplet are analogous to the case 
for an s to p transition in a spherically symmetric atom: 
light polarized along the 2 axis (tt^- light) will excite the 


state To, while circularly polarized light x±iy (cr±-light) 
excites the states T±, respectively. We show the full se¬ 
lection rules, including the Clebsch-Gordon coefficients 
in Fig. [2] 

To calculate the dipole matrix elements we use the re¬ 
cently developed variational solutions for the Is states 
for each donor from Ref. [ 18 j 3 ^ These variational solu¬ 
tions, together with the approximate form for U cc , give 
the correct energies of the six Is states for the four 
group-V donors and give good agreement with the mea¬ 
sured hyperfine coupling of the IsAi ground stated The 
results are shown in Table [T] for k l isA 1 ,isT 2 defined by 
Aab = § I2i=x,y,z l(%'l a )] 2 - Tllis quantity is related to 
the oscillator strength for the a —>■ b transition 

2t77,* 

fab=^(E b -E a )fi 2 ab , (15) 


where m* = 3 + 2/m±)~ 1 is the average effec¬ 

tive mass for the Si conduction band. In Table |T] we 
tabulate p, a b using measured oscillator strengths and 
Eq. (15). The oscillator strengths were taken from pre¬ 
viously reported values based on absorption measure¬ 
ments in doped samples!^! We also show the theoreti¬ 
cal and experimental values for the Is to 2p 0 ,± transi¬ 
tions, which are about three orders of magnitude larger. 
The forbidden IsMi to 1sT 2 transition has only been di¬ 
rectly o bserve d in absorption measurements on Bi doped 
samples which is consistent with the expectation 
that Bi has the largest overlap of the Is states with the 
nucleus. (Note that this in contrast to the deep chalco- 
gen donors in silicon, where the IsAi to 1sT 2 transitions 
are more readily observablePS) However, the measured 
dipole moment for Bi donors is 20 times larger than what 
we calculate from the variational wavefunctions. It is 
possible that Umklapp processes, neglected in the EMT, 
strongly contribute to this transition or that the varia¬ 
tional wavefunctions are inaccurate. With all this taken 
into account, we conclude that, due to the combination of 
a large spin-orbit splitting and dipole element, Bi and Sb 
are the most promising donors for the purposes of direct 
optical spin manipulations through the 1sT 2 state. 
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Finally, we end this section by remarking that for all 
the donors it would be possible to resonantly enhance 
this transition using strain and electric fields via the tech¬ 
niques described in Sec. |IVB| This would allow resonant 
manipulation via a Raman transition from IsAi to a hy¬ 
bridized state of 2p 0 and 2s and then to the 1sT 2 states. 
As the method is similar to what we describe below, we 
leave the detailed analysis to Sec. |IVB| 


(a) J = 3/2 mj = 

J — 1/2 mj = 1/2 



14 .) It) 


3/2 (b) High Field |t) 

ui 0 > AI Rfi 


(c) 

t 

I 


Low Field 

<v 0 <K AI 

T 'Ll 

mj 


mj = 

:*Ru 

= I 
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B. Electron-Nuclear Spin Polarization 


The selection rules shown in Fig. fallow optical pump¬ 
ing of any donor into an electron spin-polarized state by 
driving the system with circularly polarized light. This 
process is illustrated in Fig. [3ja— b) for x + iy polarized 
lightP^ In this section, we analyze this process in more 
detail and find the pumping rates and total electron- 
nuclear spin polarization as a function of the spin-orbit 
coupling. We consider the high and low magnetic field 
regimes, which are defined by the condition that the elec¬ 
tron Zeeman energy is much larger or smaller than the 
hyperfine splitting, respectively. 

In the high field regime, the electron and nuclear spins 
decouple and we can describe the electron spin dynamics 
independently of the nuclear spins. We take the control 
field on resonance with the J = 1/2 manifold. Assuming 
the electron initially starts in the down state, for Rabi 
frequencies fi = yiE/h -C 7 (E is the electric field ampli¬ 
tude) and short times, we can make the approximation 


W « U) + c — 11/2,1/2) + c+13/2,1/2), (16) 

c + = — (7 + i\/h)c+ + iSl/y/3, (17) 

c_ = —7 c_ + iy/2/3£l, (18) 


where the eigenstates | J, to) are given in Eq. 11 and c± 
fl /7 <C 1. Then the quasi-steady state is given by 


1 \ 

I) H- 

7 


B7 + iA/ V + i)-^A/-|r„t> 


fry + iX 


3 try + i\' 


(19) 

The optical pumping rate out of the state | j,) into the 
state | f) is 


R n = 2 7 |( T 0 t |^)| 2 


4 0 2 


A 2 


9 7 (Try ) 2 + A 2 


( 20 ) 


Therefore the minimal requirement to have rapid optical 
pumping is that A > ftry ~ (0.01 — 0.1) meV, which is 
satisfied for all the donors as shown in Table [I] 

If we allow for electron spin relaxation at rate l/Tf <C 
R jq, then this optical pumping process leads to an elec¬ 
tron spin polarization given by 


P = 


1 + Ry\Ti 

2 + Ry\Ti 


1 - 


-1 


( 21 ) 


For Bi with resonant, far-IR light at an intensity of 
10 W/cm~ 2 , y, = 1 DpS and 7 / 2 - 7 / = 1.5 GHz, we find 


FIG. 3: (a) Electron spin polarization due to optical pump¬ 
ing with cr_|_ light. The spin-flipping optical pumping term 
(dashed line) requires the spin-orbit coupling to be larger than 
the 1sT 2 linewidth. (b) In the high-field regime the electron 
and nuclear spin are decoupled. The electron spin is polarized 
into the up state (red) at the rate [see Eq. (201]. (c) I11 the 
low field regime the electron and nuclear spins are admixed 
by the hyperfine coupling (vertical arrows). This competes 
with the optical pumping to drive the system into the fully 
polarized electron-nuclear spin state (red). 


n/2ir = 200 MHz and R-yi/2n = 80 MHz. The T\ depo¬ 
larization time varies from ss 1 s in nanodevices at large 
magnetic fields to > 1000 s in bulk samples,®^®, both 
much slower than the polarization rate. Consequently, 
this optical pumping process would lead to rapid polar¬ 
ization of the electron spin p > 0.999 in ~ 100 ns. For P, 
As, and Sb, we suspect similar dynamics to Bi because 
our estimate for the dipole moment is the same order of 
magnitude (see Table []]); however, unlike Bi, these tran¬ 
sitions have never been directly observed in these donors. 

In the low magnetic field regime, we also need to take 
into account hyperfine coupling. The Hamiltonian for 
the electron-nuclear spin system is of the form 

Hhf = HujqSz — Hu} n I z + HA s ■ /, (22) 

where s(I) is the donor electron(nuclear) spin operator, 
I = 1/2, 3/2, 5/2, and 9/2 for P, As, Sb, and Bi, respec¬ 
tively, Wo( ra ) is the Zeeman energy of the electron(nuclear) 
spin, and A = (117.53,198.35,186.80,1475.4) MHz are 
the hyperfine coupling constants for (P,As,Sb,Bi). The 
two fully polarized states | fj-Oi I — -0 are eigenstates 
for any magnetic field. 

Taking the same configuration for the optical driving 
field as the high field regime we can assume Ry± <C A. In 
this case, the nuclear spin will rapidly polarize with the 
electron spin because the only pure electron spin up state 
is | 1 1) ■ The other to./ states are admixed with the spin- 
down state through the hyperfine coupling and are not 
steady states [see Fig. [3])b)] . In this limit, the occupation 
probability p m of the states | m) is perturbatively sup¬ 
pressed in l/i? n Ti, i.e., p m ~ Pm-iRn^- Therefore, 
the electron-nuclear spin polarization pi is similar to the 
case without the nuclear spins, and scales as 

( 23 ) 

2 + Kyili 

Similar to the high field regime, this allows for rapid 
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polarization of the donor electron-nuclear spin system 
pi > 0.999 in > 100 ns. 

We end this section by noting that strain will shift the 
relative energies of the six valley states, which will admix 
the states 13/2, ±1/2) with 11/2, ±1/2). When this strain 
coupling is much larger than A, then the polarization 
process will no longer be effective because the eigenstates 
become pure spin states. Using the values for the strain 
parameters in Si from Ref. |46j we find that for P and As 
the strain coupling is equal to A for strains around 10 -5 , 
while for Bi, this occurs around 10~ 4 . Strains as large 
as ~ 10 -3 are common in Si nanodevices^ 7 indicating 
that this polarization process is most applicable to Bi 
in nanodevices, but will also be achievable in low strain 
environments, such as bulk samples, for the other donors. 


IV. NEAR-IR SPIN-VALLEY-ORBITAL 
CONTROL 

In this section we consider near-IR control of the 
donors using two-photon transitions through the spin- 
3/2 donor bound exciton states D°X (see Fig. [I]). We 
calculate the dipole matrix elements for transitions to 
even-parity states and consider schemes using external 
electric fields and strain to couple to odd-parity states of 
the donor. The D°X state forms when exciton binds to 
the neutral donor state D°. Despite the indirect band- 
gap of Si, this transition is optically active because the 
point-like nature of the donor enables momentum con¬ 
servation during photon absorption or emission.^ 
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As 

Sb 

Bi 

Ed°X ~ ElsAi ( e V) a 

1.150 

1.149 

1.150 

1.147 

T (/xs)“ 

0.272 

0.183 

- 

0.0086 

f L lsA 1 ,D°X (D) 

0.033 

0.039 

0.033 

0.058 

TV (ps) C 

fo/fa 

1100 

750 

1100 

350 

Exp. 

0.24 d 

0.14 d 

0.32 d 

- 

Theory 

fo/fa 

0.06 e 

" 

" 

“ 

Theory 

/V/A 

0.03 e 

“ 

" 

" 

Theory 

0.02 e 

- 

- 

- 


a Exp., Ref. [2] 

^Exp., Ref. l48l 
c Estimated from fii s A 1 o°x- 

^Estimated from the ratio of photoluminescence intensity for no¬ 
phonon (P, As) or transverse-optical phonon sideband (Sb) transi¬ 
tions ending in Is versus 2s reported in Ref. 1491 
"Theory, this work (neglects central cell corrections). 

TABLE II: Table of physical parameters relevant for near-IR 
spin-valley-orbital control. As the free exciton recombination 
energy is around 1.1545 eVp^we can see that the D°X state 
has a binding energy around 5 meV (note the Si conduction 
band gap E g fa 1.17 eV), r is the D°X state lifetime limited 
by Auger scattering, and r r is the radiative lifetime calculated 
from fi la A 1 ,D°x- 


In Table [TT] we compile some of the relevant parame¬ 
ters for near-IR control of the donors. It is important to 
note that, due to Auger recombination processes, these 
states are far from radiatively broadened and the ratio 
of their natural linewidth to the radiative linewidth is 
r r /r ss 4000 for P and As and 40 000 for Bi. As a re¬ 
sult, creating an efficient optical interface to the donor 
spin states with these states is challenging. They have 
proven to be a powerful resource for electron-nuclear spin 
to charge conversion for spin readout, spin control, and 
state irhtializationP ^ 12 f l4 l Here we explore their potential 
use for quantum control of the spin-valley-orbital states 
of the donor. 

To characterize these optical transitions it is impor¬ 
tant to consider how the Coulomb interactions between 
electrons and holes affect the D°X ground state. To treat 
this problem we use a Hartree-Fock approximation with 
the assumption the two electron ground state is predom¬ 
inantly a spin- singlet in a single orbital state with A\ 
symmetryp 4 50 This is a good approximation because of 
the large orbital splittings of the donor. For valley elec¬ 
tron and hole envelope functions F and Fh, respectively, 
the Hartree potentials take the fornPdL- 0 - 


U e H 

UhH 

V e (r) 

V h (r) 


U{r) - eV e (r) - eV h (r), 

4^ + 2el4(r) . 


e 1 


47T6oe 6 
e 


X 




\r — r' 


47reoe 




(24) 

(25) 

(26) 
(27) 


Within the single-valley EMT, this gives rise to the 
Hartree-Fock (HF) equations for the two electrons and 
the hole 


£ F H ~ 

C e ± ev 

(■t + U eH )F%, 

(28) 

z h Fh = 


(29) 


where the index v refers to the valley state and rrih = 
0.49(0.16)m 0 is the heavy (light) hole mass. 

To solve the HF equations we use a pair of variational 
solutions for the electrons and the holes. For the elec¬ 
trons we take a Kohn-Luttinger form for the variational 
wavefunction shown in Eq. ([8]) , while for the hole we take 
a hydrogenic wavefunction of the fornP^l 


F, 


H 


_i_ r —e~r/2c 

2 C 5/2 ^3/ 


(30) 


with c a variational parameter. F^ vanishes at the ori¬ 
gin, where the donor potential is repulsive, and is spher¬ 
ically symmetric, which is consistent with the symmetry 
of the valence band and the symmetric form of UhH- To 
find the variational parameters we fix c and minimize 
the expectation value of Eq. (28) with respect to a and 
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FIG. 4: Selection rules, Clebsch-Gordon coefficients, and overlap factors for exciting from the IsAi and 2sAi states to the D°X 
state for (a) left-hand-circularly polarized light cr_, linearly polarized light along the z-axis n z , and (c) right-hand-circularly 
polarized light cr+. 


b , and similarly for the hole, we fix a and b and mini¬ 
mize the expectation value of Eq. (29) with respect to 

c. Self-consistent variational solutions occur when these 
two minimization procedures produce the same values 
for a, b , and c. For the heavy hole we find a = 2.79 mn, 
b = 1.60 nm, and c = 1.34 nm, while for the light hole 
we find a = 2.74 nm, b = 1.57 nm, and c = 1.99 nm. 
This is to be compared to the case for the D° state 
where the Kohn-Luttinger values are a = 2.51 nm and 
b = 1.44 nm.® 


A. Two-Photon Transitions to Even-Parity States 


We now explore the potential for spin-selective Raman 
transitions from the \sA\ ground state to the even parity 
excited orbital states of the donor such as the 2sAi and 
3sA! states. The two electrons in the D°X ground state 
are expected to primarily pair in a spin-singlet in the 
\sA\ orbital, as a result the nsAi orbitals will become 
admixed with the ground state through the Coulomb in¬ 
teractions between the electrons and holes. Such two- 
photon transitions have been used to explain some of 
the satellite emission lines observed in photoluminescence 
spectra of the donor bound excitons.’ 2 ® In Table |Il] we 
list the ratio of the dipole moments for the no-phonon 
transitions D°X —> lsA\ vs the D°X —> 2sAi states 
based on these early measurements.® 

We can also calculate this ratio using the Hartree-Fock 
theory. Our variational solution has the same symmetry 
as the ns states of the singly occupied donor, which allows 
the expansion 

IV’S) =X)/ 3 "i ( 31 ) 

n 

Pn = (VnsliO = J drF"* s {r)F»{r). (32) 


We approximate the ns states by an orthogonal set of 


hydrogen-like wavefunctioni 51 -^ 

Ks = r-KnrX U 2p/n)e-r/\ (33) 

V 7r n°a z b 

p = \J (x 2 + y 2 )/a 2 + z 2 /b 2 , (34) 

where L^}_ 1 are the generalized Laguerre polynomials for 
the ns states and we take the same values for a and b 
as the Is state to ensure orthogonality. The ratio of the 
dipole matrix elements is then given by 

Pn'sAi,D°X Pn' (35) 

PnsAi,D°X Pn 

In Table El we tabulate these ratios for the first few Is 
states. In general, we find /3 n ~ 1 jn as n increases. 
We show the full selection rules, including the Clebsch- 
Gordon coefficients, for these transitions in Fig. [4] 
Additional two-electron transitions have also been ob¬ 
served to the 1 sE and 1sT 2 stat esP However, these states 
become admixed with the D°X ground state through 
valley-orbit interactions, which we have not included in 
the Hartree-Fock analysis. 


B. Two-Photon Transitions to Odd-Parity States 


In this section, we show how to achieve two-photon 
transitions to the np-states by applying either a static 
electric field or strain to the donor system. This addi¬ 
tional control is required because the hydrogenic model 
is a good approximation for the donors and parity is a 
good quantum number. Parity can be broken by apply¬ 
ing a static electric held, which will mix ns and np-like 
states. We analyze this case in Sec. |IV B 1| below. Al¬ 
ternatively, valley-orbit interactions can break the parity 
symmetry because the tetrahedral group does not con¬ 


serve parity. Similar to the case in Sec. Ill there is a 
direct transition from the D°X state to the 2p° states 
which are antisymmetric combinations of the same val¬ 
ley states. This occurs with the relative dipole moment 


P2p 0 ,D°X 

PlsA 1 ,D°X 


i / dPru*_ z u z F2p F^ z s\nk 0 z, (36) 





































Unfortunately, this is suppressed by a factor of ~ 10 3 
compared to the nsA\ transitions, due to the rapidly 
varying Bloch phase. In Sec. |IV B 2[ we show that these 
transitions are resonantly enhanced for certain values of 
strain when the 2 s states are degenerate with the 2 p 
states to within ~ (10 — 100) /xeV. 

1. Static Electric Field 

When an electric field E is applied to the donor the 
single-valley EMT is modified to 



EF^r) = + U(r) + eE ■ r F M (r). (37) 


By analogy with Stark effect for a hydrogen atom, we ca n 
write an ansatz for the 2po state in the ±fco 2 vallejIEU 

= Nz{ 1 + g: ) e -V( i WSW ] ( 38 ) 

where N is a normalization constant, q 7 a, and b are vari¬ 
ational parameters and we took E parallel to the z -valley 
axis. For each value of |E we minimize the energy expec¬ 
tation value to find the variational parameters. For finite 
q the symmetric combination of these valley states has a 
direct overlap with \sA\ Hartree-Fock orbital, which re¬ 
sults in the ratio of dipole moments 

LL2 Po z,d°x = J_f d 3 rF~;F ( 39 ) 

UlsAuDOX V3 J P ° 

For small fields, we are justified in taking the zero-field 
envelope function for the Hartree-Fock orbital because of 
the weak dipole moments of the D°X state compared to 
the 2p 0 state. 

In Fig. [5] we show this ratio under the application 
of experimentally relevant electric fields (note the the 
donors are ionized at electric fields around 2 U//nrP^). 
We also plot the relative shift in the energy of the 
2p 0 state and the D°X state with applied electric held. 
For the D°X we used the Stark shift parameter 2 p$ ss 
33 peV/(V/p m) _1 measured in Ref. [14] The shift in the 
D°X state is just a few-percent of the binding energy, in¬ 
dicating the zero-held envelope is a good approximation 
to the wavefunction. The shift in the 2po energy is much 
larger, but it is still small compared to the 2po binding 
energy. From the figure we see that at the optimal value 
of the electric held, the ratio of the dipole moment to the 
odd-parity 2 p 0 state is the same order as the ratios we 
found for the even-parity transitions at zero-held. Simi¬ 
lar arguments also hold for the various 2p± states. Thus 
we can conclude that a static electric held is a realistic 
approach for achieving two-photon transitions to these 
odd-parity valley-orbit states. 


FIG. 5: (Black) Relative dipole moment for optical transitions 
from D°X ending in the odd-parity 2po state with increasing 
electric held, which breaks the parity symmetry. (Blue/Red) 
Energy shifts of 2po/D°X states relative to their binding en¬ 
ergy with applied electric held. The peak in the dipole mo¬ 
ment occurs when both shifts are relatively small. For small 
helds the dipole moment is small because parity is only weakly 
broken, while for large electric helds the 2po state becomes 
further shifted away from the donor site where it has a large 
overlap with the lsA\ state. 


an additional electric field. In Fig. [6] we show the shift 
in the energy levels of the D° and D°X states with a 
compressive pressure along the [001] axis for phospho¬ 
rous donors using the model developed in Ref. [HI For 
pressures around 5(30) MPa we see that the lower(upper) 
valleys of the 2s states becomes resonant with the up¬ 
per (lower) valleys of the 2Po(±) states. In Fig. [61(b) 
we also show the 2 sAi state for bismuth donors, winch 
has the coincidence that these resonances occur near the 
same value of the pressure around (15-20) MPa. 

Within the single-valley EMT these states remain de¬ 
coupled even under the application of a small electric 
field because the electron is localized in opposite valleys. 
However, when the valley-orbit interactions are included, 
these become avoided crossings because the donor po¬ 
tential can induce transitions between valleys. Looking 
more closely at the 2 s l — 2p% crossing (£(u) refers to the 
lower(upper) valley states under the application of strain) 
we can express these states as 


a 2s e ,± 


a 2 Pq,±x 


a 2p%,±y 


-^(0,0,0,0,1,±1), 

(40) 

^=( 1 , ± 1 , 0 , 0 , 0 , 0 ), 

(41) 

^=(0,0,1, ±1,0,0), 

(42) 


The valley-orbit interactions are given by the matrix el¬ 
ements 


2. Strain Enhanced Valley-Orbit Interactions 


= {ip n \U\i/)' n ). (43) 


We now show that it is possible to achieve comparable 
two-photon transitions to the 2 p states without applying 


Because of the odd-parity of the 2 p 0 states and the even 
parity of the 2s states, this will only be nonzero for the 
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FIG. 6: (a) Spectrum of donor transitions for Si:P from the 
D° ground state to the excited D° states and the D°X state 
with increasing compressive force along the [001] axis. Note 
for P the spin-orbit splitting of the 1 sT 2 and 2 p states is neg¬ 
ligible compared to the linewidth of these states, (b) Ener¬ 
gies of 2 p and 2s states with energy measured relative to the 
conduction band at zero detuning. Resonances appear at spe¬ 
cific strain values. Also shown are the states for Si:Bi where 
resonances appear between 2s and 2po and 2 p± at similar 
pressures, (c) Avoided crossing between the symmetric com¬ 
bination of lower-valley 2s l states and the four upper-valley 
2po states for A = 10 /reV. Two of the five states repel each 
other due to valley-orbit interactions, while the other three 
are unaffected. 


states a 2s ‘,+ and a 2p ^- x (- y ), in which case 


A = 2 i d 3 r ulu x F 2 *(r)U(r)F 2 po (r) cos k 0 z sin k 0 x 

(44) 

A simple estimate suggests this should be ~ 10 3 times 
the valley-orbit splitting of the Is states, which implies 
that it can be as large as (10 — 40) /zeV, depending on 
the donor. This is comparable to the lifetime of these 
states. In Fig. [HJc) we show a small region of pressure 
around this avoided crossing for phosphorous. We see 
that the width of the avoided crossing corresponds to 
about 1 MPa, indicating that it would be possible to 
stabilize the system at this point with suitable control of 
the stress. At the avoided crossing when E 2s e <+ = E 2p u 
the eigenstates are given by 


1 


AS = A :a =-^(1,-1,1,-1; 1,1) 


(45) 

A£ = 0:a = {-^(1,1,0,0;0,0), (46) 

-^=( 0 , 0 , 1 , 1 ; 0 , 0 ), ^( 1 ,- 1 , - 1 , 1 ; 0 , 0 )}, 


AE=~A: a =—(—1,1,-1,1; 1,1), 


(47) 


where A E = E — E 2p <* and it is should be understood 


that the first four entries are associated with the 2pg val¬ 
ley envelopes and the last two entries with the 2s l valley 
envelopes. At this resonance, the ratio of the dipole mo¬ 
ments will then be given by 


F2p 0 ,nox _ 1 f 

„ ~ o / a 1 r 2s r ez ~ o o 1 

PlsA 1 ,D°X 3 J O pi 


(48) 


where we took the zero-stress form for the Hartree-Fock 
solution, which is valid approximation for these small 
stresses because the ground state wavefunction is pro¬ 
tected from deforming by the large valley-orbit splitting. 
This results in an enhancement of the dipo le moment by 
a factor of ~ 10 3 compared to the Eq. (36). 

We end this section by noting that although strain 
is common in silicon nanodevices (it can be as large as 
~ 10 -3 corresponding to P w 200 MPeP^), the approach 
described here would require additional static tuning of 
the strain to bring the donor states near these resonances. 
Furthermore, in this discussion we have only considered 
the influence of axial strain, in realistic nanodevices there 
is also a contribution from shear strain. Shear strain 
leads to additional energy shifts and it changes the effec¬ 
tive masses of the conduction bandP^ This latter effect 
will change the valley splittings, but will not qualitatively 
change the character of these resonances. Therefore, the 
main challenge for exploiting these resonances in nanode¬ 
vices will be achieving the resonant strain condition in a 
deterministic manner. Alternatively, mixing between s 
and p states might be induced by other effects such as 
proximity to an interface. In bulk samples, these reso¬ 
nances should be more readily achievable by applying a 
global stress to the sample, as we have assumed in this 
section. 


V. CONCLUSIONS 

We have shown how to use far-IR or near-IR optical 
fields to achieve spin selective excitation of the group-V 
donors from the ground state to excited orbital states. 
In the case of far-IR light, we calculated the dipole mo¬ 
ments for transitions which are parity-forbidden within 
the hydrogenic approximation for the donors, but are al¬ 
lowed due to valley-orbit interactions. These transitions 
have only been directly observed in the case of the Bi 
donor. For Bi, our calculated dipole moment disagrees is 
two orders of magnitude smaller than the estimate based 
on these absorption measurements. This suggests that 
further experimental and theoretical work is required to 
account for such a large discrepancy. We then showed 
that these transitions obey selection rules consistent with 
excitation from a J = 1/2 ground state to spin-orbit split 
J = 3/2 and J = 1/2 excited states. This enables op¬ 
tical pumping into electron-nuclear spin polarized states 
through the application of circularly polarized light. We 
calculated the timescale and final polarization for this 
process and found it is most promising for Bi donors, due 
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to their large spin-orbit splitting and strong valley-orbit 
interactions. 

For near-IR light one can make use of the J = 
3/2 donor bound exciton state D°X, which has weak, 
dipole allowed transitions from the donor ground state 
to D°X. We used a Hartree-Fock variational solution to 
the electron-hole wavefunctions for this state. Using this 
solution we estimated the relative strength of the transi¬ 
tion from the lsA\ ground state to D°X versus the tran¬ 
sitions from the nsAi excited states to D°X. For the Is 
versus 2s transitions, our calculated values agree with 
previous experimental measurements within a factor of 
five. Resonant excitation on these two transitions would 
enable spin-selective, two-photon, A-transitions from the 
ground state to the nsAi excited states. We then showed 
that one can achieve similar A-transitions to 2 p states by 
breaking the parity symmetry present within the hydro- 
genic approximation for the donors. We considered two 
approaches. In the first approach, an applied electric 
field directly breaks the parity symmetry and admixes 
the nsA-[ states with 2 p states. The second approach uses 
strain to bring the 2s and 2 p states into resonance. In this 
case, the valley-orbit interaction is resonantly enhanced 


and the antisymmetric combination of 2 p valley states 
becomes strongly admixed with the symmetric combina¬ 
tion of 2s valley states. This enhances the A-transition 
from the IsAi to the 2 p states by a factor of ~ 10 3 at 
zero applied electric field. 

We anticipate that these results will help guide future 
spectroscopic investigation of the donors. In particular, 
many of the effects investigated here probe the poorly un¬ 
derstood interaction between the donor bound electrons 
and holes and the donor nucleus. For quantum informa¬ 
tion applications, these results provide a path forward for 
achieving full optical quantum control and state prepa¬ 
ration of group-V donors, as well as a coherent interface 
between donor spin qubits and single photons.^ 
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